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Abstract 

Many multi-dimensional consistent discrete systems have soliton solutions with nonzero back¬ 
grounds, which brings difhculty in the investigation of integrable characteristics. In this letter we 
derive infinitely many conserved quantities for the lattice potential Korteweg-de Vries equation. 
The derivation is based on the fact that the scattering data a(z) is independent of discrete space 
and time and the analytic property of dost solutions of the discrete Schrodinger spectral prob¬ 
lem. The obtained conserved densities are different from those in the known literatures. They are 
asymptotic to zero when \n\ (or \m\) tends to infinity. To obtain these results, we reconstruct a 
discrete Riccati equation by using a conformal map which transforms the upper complex plane to 
the inside of unit circle. Series solution to the Riccati equation is constructed based on the analytic 
and asymptotic properties of dost solutions. 
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1 Introduction 


For a nonlinear partial differential equation, possessing infinitely many conservation laws indicates the 
equation is integrable. In the fully discrete case, multi-dimensional consistency [l]-[^ provides a kind 
of strong integrability, under which and extra conditions lattice equations defined on quadrilateral 
were classified and discrete Boussinesq(DBSQ)-type equations were found as multi-component 
models. For these multi-dimensional consistent equations, infinitely many conservation laws were 
already constructed in several different ways [5-14 . However, these known conservation laws are only 
formal ones — the nonzero boundary condition at infinity (asymptotic behavior) of the solutions of 
the considered lattice equations was neglected. 

How does asymptotic behavior play roles? A discrete conservation law is defined in the form 


{Em ~ ^)En,m — {En ~ 1) Jn 


( 1 . 1 ) 
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where Em and are shift operators defined as Emfn,m = /n,m+i and Enfn,m = fn+i,m- We request 
Fn,m and Jn,m tend to zero fast enough when \n\ —)■ oo so that = 0 and the 

"VOO 


sumation ^n,m makes sense. As a result, ^n,m is indenpendent of m and represents 

a conserved quantity with respect to m. 

To explain more, as an example, let us look at the lattice potential Korteweg-de Vries (IpKdV) 
15,16 , also known as HI equation in the Adler-Bobenko-Suris (ABS) list [^, 


equation 


{u — u){u — u) = 


( 1 . 2 ) 


where u = Un,m: u = Un+i,m, u = Un,m+i, u = Un+i,m+i, and lattice parameters p for direction n and 
q for direction m are arbitrary constants. Unlike most of continuous integrable systems, u = 0 is not 
a solution to the IpKdV equation (1.2). It has a simplest solution 


u = pn + qm + c 

with a constant c. The V-soliton solution has the form 


u = pn + qm + c + w 


where 


w 


0 , 


n 


oo (or 


m 


oo). 


(1.3) 


(1.4) 


(1.5) 


Solution (1.3) is usually called a background of solitons. It has been known that most of members in 


the ABS list and DBSQ-type equations have either linear or exponential backgrounds (cf. [17| - [^ ). 
The nonzero background will bring difficulties in Ending convergent conserved quantities. For the 
IpKdV equation (1.2), due to the background H, the equation itself is not in a conserved form and 
the discrete integral of function of u, for example, does not make sense because with 


the background (1.3) the summation is not convergent at all. With regard to discrete conservation 
et us 

[12]), which is in the form (1.1) with 


laws, as an example, let us take a look at the second conservation law of the IpKdV equation derived 
from its Lax pair (cf. 


Er},,rn. — 


1 


Jn.m — 


1 


{u — u){u — u) 


{u — u){u — u) 


( 1 . 6 ) 


Obviously, for u defined in (1.4) with the background (1.3), neither Fn,m nor Jn,m is asymptotic to 
zero when |n| —)• 0 or |m| —)• 0. Although one can artificially add constants —l/4p^ and l/2p{p — q) 
to the above Fn^m and Jn,m respectively so that they get zero asymptotic behaviors, these constants 
can not be derived naturally. This is the drawback. 


In this letter, for the IpKdV equation (1.2), we derive its infinitely many conserved quantities of 
which each corresponding conserved density F^^m 0 when |n| —>■ oo and the conserved quantity 
X)^_oo ^n,m is convergent. Our approach is a by-product of the inverse scattering transform (1ST). 
In the 1ST, the scattering data a(A) (connected with the transmission coefficient T by T = l/a(A)) 
is time-independent and it is asymptotically related to wave functions (Jost solutions) in a certain 
way. a(A) can be expanded in the complex plane and coefficient of each term provides a conserved 
quantity (functional). Such an approach was first proposed by Zakharov and Shabat in 24 where 
they derived conserved quantities for the nonlinear Schrodinger equation. Later, it was applied to 
the Ablowitz-Ladik hierarchy [^. The approach was also used to obtain conserved quantities for 
the KdV equation (see [^). However, for the fully discrete multi-dimensionally consistent IpKdV 
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equation (1.2), particularly, u having non-zero background, we need to reconstruct a new discrete 


Riccati equation by using a conformal map which transforms the upper complex plane to the inside of 
unit circle. Series solution to the Riccati equation is constructed based on the analytic and asymptotic 
properties of Jost solutions. 

The letter is organized as follows. In Secwe recall some necessary results in the direct scattering 
problem of the IpKdV equation given in (^ . Then in Secwe derive infinitely many conserved 
quantities for the IpKdV equation. Finally, Secj^is for conclusions. Besides, we will revisit solutions 
to the Riccati equation in Appendix. 


2 Preliminary 


In this section we recall some results in the 1ST procedure of the IpKdV equation described in 22 


These results are necessary in our approach. Throughout the paper we suppose the spacing parameters 
p > 0 and q > 0. 


The IpKdV equation (1.2) has the following Lax pair 


(h 

4>2 

4>2 


—u UU+ p^ — 
— 1 u 


—u UU + 


-1 


u 


(t>i 

4>2 

4>i 

(f>2 


One can rewrite (2.1a) into a scalar discrete Schrodinger problem |27[|28| 

g- sg + {p^ + z^)g = 0 , 

where 


S — Sn.m. — U 11 


(2.1a) 

(2.1b) 

( 2 . 2 ) 

(2.3) 


and we have replaced the wave function 4)2 with g and the spectral parameter r with iz. It then follows 

(2.4) 


from the structure (1.4) that 


s ~ 2p, |re| —>■ oo. 

This is the boundary condition for s or tt, and we request s falls in the space (cf. [22|[^) 


P,= 


Sn,m ■ ^ ^ 1^^ (Sn,m 2p) | < OO^ , fJ, — 0,1 

n=—oo 

We note that in direct scattering problem m is treated as a constant or zero. 


2 . 


(2.5) 


Under the boundary condition (2.5), the difference spectral problem (2.2) admits Jost solutions ip 
and with asymptotic property 


ip{n] z) ^ {p — iz)"', (n —>■ —oo), 
z) ^ {p + iz)", (n—)-oo). 


(2.6a) 

(2.6b) 


where p{n; z) and ^(n; z) are analytic in the upper complex plane Imz > 0 and continuous in Imz > 0. 
Besides, 


ip{n;z){p-iz) "l + o(^^y (lkl|-^oo). 


(2.7a) 
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(2.7b) 


'tp{n; z){p + iz) + (|| 2 ;|| —> cx)). 

There is a linear relation 


ip{n\ z) = a{z)'il){n-, —z) + b{z)'il^{n; z), {z G M). 


a{z) is independent of n and expressed as 

a{z) = ^Wnipin] z),'ip{n; z)], 

where the Casorati determinant with respect to n is defined as 


Wn[/(n),5r(n)] 


+ z^)-^ 


f{n;z) 
f{n + l-,z) 


g{n;z) 
g{n + l-,z) 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


Thanks to the analytic property of ip and ip, the relation ( |2.9| ) can be extended to the upper complex 
plane, i.e. a{z) is analytic in Imz > 0. For the asymptotic property, a(z) ~ 1 when ||z|| —)■ oo and 
a{z) can be expanded in the neighbourhood of oo as 




i=i 


( 2 . 11 ) 


a{z) has finite zeros {z = and due to analytic property and uniqueness, a(z) can be 

expressed as (cf. [^) 


N 

/ \ TT ^ , 

=n 

j=l 3 

Therefore a{z) is independent of both n and m. 


3 Infinitely many conserved quantities 

3.1 General scheme 

Making use of the relation (2.9) and the asymptotic behavior (|2.6|), we hav^ 


az ~ 


p{n]z) 


p + iz — 


ip{n + 1 ; z) 


and consequently, 
where 


2iz{p — iz)^\ ip{n] z 

Ina{z) ^—ln2iz + Gn + Bn, (n —)> oo), 
p{n-,z] 


, (n->oo). 


Gn = In 


{p — izy 


Bn = ln(p + iz - 9) 


*In the case of the continuous KdV equation, the counterpart of (13.11) is 




However, the term on the r.h.s was missed in 


26 


(3.1) 

(3.2) 

(3.3) 
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and 


9 = e{n) = 


(p{n + 1; z) 
ip{n;z) 


Noting that lim„_s>_oo = ln2iz, we have 

hia(z) ~ Gn +-Bn --B_oo, (n^oo). 

Here B-oo stands for hmn->.-oo Bn and hereafter we also employ similar notations such as 
G-oo without confusion. Next, defining 


Fn = (Bn - l)Gn = In 


9 


p — IZ 


and noting that G-oo = 0 due to the asymptotic property (2.6a), we hnd 

n—1 n—1 

lna( 2 ;)~ (n-^oo). 

j=—oo j=—oo 

Thus, since a(z) is independent of n it is then recovered via 

n—1 n—1 


hia(z) = lim 

n^OQ 


[ ^ ^ (Bj+i Bj 

j=—oo j=—oo 


9 


— ^ ' In . + Boo B—oo- 

p — iz 

n=—oo 

In principle, this will provide a set of conserved quantities {Kj'\ by expanding 


Ina(z) = In^l + ^ 
i=i 

3.2 Riccati equation and conformal map 


oo oo 

a 


E 




z^ 


The known discrete Riccati equation w.r.t. 9 corresponding to (2.2) (with g = (p{n]z)) is 

99 — s9 + + z^. 

It admits a series solution |12l 

O OO 


9 OO 




with 


in which 


1 ^ ~ 

^i+i = j = 0,l,2,---, (00 = 1), 


i=0 


01 = 4 , 02 = 4 ( 4+4 

SS SS \SS gj 


(3.4) 

(3.5) 
Boo and 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11a) 

(3.11b) 

(3.12) 
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Obviously, with the above form (3.11) for 6, it is hard to expand In into a series of z with 


negative powers which matches the expansion (3.9) 
To over come this obstacle, we consider 


9 


UJ = 


p — IZ 


under which 


r„ = lnw, Bn = in(-—- Lv) +ln{p - iz), 
\p — iz / 

and the Riccati equation w.r.t. oj takes the form 

su p + iz 


UJUJ -^—h 


= 0 . 


p — iz p — IZ 


To solve it, we introduce a conformal map 


A = 


p + iz 
p — iz 


(3.13) 

(3.14) 

(3.15) 

(3.16) 


to transform the upper z-plane to the inside of unit circle on the A-plane, as depicted in Figj^ 



Figure 1: Conformal map (3.16) 
Under this conformal map and noting that 


z = ip 


1-A 

TTa’ 


(3.14) reads 


r„ = Inw, Bn = ln(A — a;) + In 


2p 

TTa’ 


(3.8) reads 


1 »\ 

lna(ip ^ + a ) ^ ^ Incj + Roo-R- 


and the Riccati equation (3.15) is written as 

u(uj — x) — XxLo + A = 0, 

where 

X = s/2p = (u — u)j2p, 


(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 
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which tends to 1 as |n| —?• oo. (3.20) has the following solution, 


UJ = X 


1 + 

j=i 


< 1 ) 


with 


wi = 1-, 

X X 

V XX/ XX 
j 


(3.22a) 

(3.22b) 

(3.22c) 

(3.22d) 


CUj-i-i — Wj ^ ^ OJk ^j+l—kj j — 2,3, • • • . 
k=l 

Noting that the fact x ~ 1 as |n| —)> oo, we immediately have the following. 

Proposition 1. All the {ujj} defined in (3.22) are asymptotic to zero as \n\ — >■ oo. 

3.3 Infinitely many conserved qnantities 

From now on everything is considered on A-plane. Based on proposition for the term Bn in (3.18) 
we find B^o = B-^o, from which and (3.19) it follows that 


In a I ip 


1 - A 
'l + A 


^ In 


ui. 


Theoretically, the l.h.s. has the expansion (at the origin point of A-plane) 

1 - “ 

i=i 


In a (ip 


1 +A 


(3.23) 


(3.24) 


due to the property that is analytic inside the unit circle, while the r.h.s. of (3.23) can be 

explicitly expressed as 

OO OO OO OO 

Inw = Inx-I- ln(l -I- yy cjjA-^ 


n=—oo 
oo 


00/00 


yz In^ + Z] Z] M'- 

t=—OO j=l \ n=—oo J 


(3.25) 


Here, {hj{u:)}, we call them /i-polynomials for short, are the polynomials defined in the following way 
(cf. [p]) 

y OO \ OO 

(3.26) 


ln( 1-P y^^WjAM =y^/ij(cj)A/ 

^ i=i '' j=i 

where to = {uji,ui 2 , • • •) and the first three of {hj{u/)} are 

hi{u) = cji, h2{uj) =-^Lof + aJ2, hfita) = - ujiuj 2 + u}3. 

A general description for such polynomials is given in Appendix 

Thus, it is clear that the infinitely many conserved quantities are the following. 


(3.27) 
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Proposition 2. The infinitely many eonserved densities of the IpKdV equation (1.2) are 


Ho= \nx, Hj= Y hjiu:), 


(3.28) 


and the eorresponding conserved densities are 

Fo = lnx, Fj = hj{u), 

where x = {u — u)l2p and the h-polynomials are defined in Appendix\^ 

For these conserved quantities or conserved densities, we have further results. 


(3.29) 


Proposition 3. Under the criteria given in lOl, the eonserved densities (3.29) are non-trivial and 
not equivalent to each other. 


This can be proved as done in by making use of /i-polynomials to calculate orders of conserved 
densities. 


At the end of the section we note that it is not necessary to compare the conserved densities (3.29) 
with 

Fq = - \n{u-u), F' = /ij(0), j = l,2, ••• (3.30) 


which were derived in 12 based on the series solution 9 given in (3.11). Here 0 = {9i,92,-' ')■ In 
fact, on one hand, Fj are functions of {un+i_s} while Fj are functions of {u„+s} where s = 0,1, 2, • • •. 
They can not be unified. On the other hand, any series solutions to the Riccati equation ( |3.10| ) should 
agree with known analytic properties. In this sense, the expansion (3.11a) is wild. 

A further discussion on the series solutions to the Riccati equation (3.20) is given in Appendix [B| 


4 Conclusions 


We have derived infinitely many conserved quantities for the IpKdV equation (1.2). Different from 


those known ones, these conserved quantities are convergent thanks to the zero asymptotic behaviors 
of conserved densities. They are also nontrivial and not equivalent to each other. The derivation 
is based on analytic and asymptotic properties of dost solutions of the discrete Schrodinger spectral 


problem (2.2). The conformal map (3.16) is employed to convert the upper z-plane to the inside 


unite circle of A-plane. Such a spectral transformation enables us to have a new series solution of the 
Riccati equation with respect to A and the solution agrees with the analytic and asymptotic properties 
of dost solution ip{n; z). This new solution leads to desired infinitely many conserved quantities, which 
presented in Proposition]^ Another set of conserved quantities with respect to (m, q) can be obtained 
by replacing {n,pf') with {m,qfi). 

It has been seen that the derivation of conserved quantities or conservation laws for discrete systems 
is more complicated than continuous case. One reason is due to the obstacle of non-zero background 
of soliton solutions in discrete case. It is specially important to construct conserved densities which 
are asymptotic to zero when solutions have non-zero background. Another point is Riccati equation 
and its solutions. For the IpKdV equation, we have three versions of its Riccati equation, which are 
(3.10), ( 3.15[ ) and ( |3.20[ ). Even for (3.20) that we finally used in the paper, it has two different series 
solutions (see Appendix but only one agrees with the analytic and asymptotic properties of dost 
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solution ip. We hope by the discussion of the letter we can emphasize the importance of asymptotic 
property and analytic property of discrete wave functions in the research of discrete systems. 

It is possible to examine asymptotic behavior of solutions of other equations in the ABS list which 
are solved via 1ST [^, and make use of analytic and asymptotic properties of dost solutions to 
derive reasonable conserved quantities. It is also possible to re-examine Lax pair approach in [12] 
and construct reasonable Riccati equations and their series solutions for those multi-dimensionally 
consistent lattice equations (see the list in [^). 
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A 


/t-polynomials 



The following expansion holds. 


2 = 1 


ln(l + ^tj/c*) = hj{t)k\ 
i=i 


where 


and 


hjit) = E (-l)'“'""(|a| - 1)!^. 


I«N=J 


t = (tl, t2, • • • )> « = (Ol, 02, • • • ); Oj G {0, 1, 2, • • • }, 

OO OO OO OO 

q! = ]^(Q;j!), \a\ = '^ai, ||q|| = 

2=1 2=1 2=1 


i=l 


The first few of {hj{t)} are 


hi{t) = ti, 

h2{t) = --tl + t2, 

^3(t) = - ht2 + h, 

hiit) = + tjt2 - hts - + ti. 


B Proper solution to the Riccati equation (3.20) 


(A.la) 

(A.lb) 

(A.lc) 

(A.ld) 

(A.2a) 

(A.2b) 

(A.2c) 

(A.2d) 


In this section, we investigate how a proper series solution of the Riccati equation (3.20) is constructed 
based on the analytic and asymptotic properties of the dost solution ip{n; z). 
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For the Riccati equation (3.20), in addition to the series solution (3.22), it has a second solution: 


j=i 


with 


^+1 - + ~ 


1 ^ 


(B.la) 


(B.lb) 


k=l 


This solution has asymptotic property w ~ A as |re| —>■ oo, but the expansion is ‘wild’ because it does 
not agree with the analytic and asymptotic properties of ^p{n;z). 

Let us explain this in the following. In the light of asymptotic behavior (2.7a) and noting that 
(|3.13|), i.e. oj = there is 

^ K’ ‘^[n]z) p—tz^ 


a; ~ 1 + O^-^, \\z\ 

which means w can be expanded at oo as 


oo. 


UJ = 


1+E 




i=i 


z3 


Meanwhile, the conformal map (3.16) indicates 


and consequently, 


1 1 

- = i(l + 2^A^'), 

z ip 

i=i 

— = 7tV(1 + 0(A)). 

z^ {ipy^ ^ 


(B.2) 

(B.3) 

(B.4) 

(B.5) 


This means, by substituting the above into (B.3), oj should have the following expansion in terms of 
A: 


and the term wq is 


OJ = 


Wo = 


Wo + WjA-^ 

(B.6) 



oo (j) 

1 + / ,. ... 

(B.7) 


This is nothing but w(n; z = ip) in light of the expression (B.3). To calculate oj{n]ip) we take z = ip 
in the spectral problem (2.2) with g = (p{n; z), we have 

ip{n + 2; ip) — 2pxip{n + l;ip) = 0, 

which gives 


/ . , ip{n + l-,ip) 1 

a;(n; ip) = . . . — = x. 

ip) Ip 


(B.8) 


Such a value, together with the expression (B.6), coincides with the expansion (3.22). 

Thus, it is clear that the expansion (3.22) is actually well-based on the analytic and asymptotic 
properties of the wave function p(n,z) while the expansion (B.l) not. In other words, the analytic 
and asymptotic properties of wave functions play important roles in the procedure to find conserved 
quantities for discrete integrable systems. 
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